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We study properties of excited states of an elongated Fermi superfluid along the BEC-BCS 
crossover including the unitarity limit. Analytic expressions for the sound velocity in an inho- 
mogeneous as well as homogeneous Fermi superfluid along the crossover are obtained on the basis 
of the hydrodynamic theory. The complete excitation spectrum of axial quasiparticles with various 
discrete radial nodes are presented. We discuss the feasibility of measuring the sound velocity and 
the multibranch Bogoliubov spectrum experimentally. 
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I. INTRODUCTION 

Strongly interacting two-component Fermi gases pro- 
vide a unique testing ground for the theories of exotic 
systems in nature. In atomic Fermi gases, tunable strong 
interactions are produced using the Feshbach resonance 
jH Q< Q- By sweeping the magnetic field in the Fes- 
hbach resonance experiments, magnitude and nature of 
the two-body interaction strength changes from repulsive 
to attractive. Across the resonance the s-wave scattering 
length a goes from large positive to large negative values. 
The fermionic system becomes molecular Bose-Einstein 
condensates (BEC) for strong repulsive interaction and 
transforms into the Bardeen-Cooper-Schrieffer (BCS) su- 
perfluid when the interaction is attractive. The first 
observations of BEC of molecules consisting of loosely 
bound fermionic atoms [a |y, |]j initiated a series of ex- 
plorations 0, H, IE OH LUj of the crossover between BEC 
and BCS superfluid. The size of fermion pair condensates 
smoothly increases from the BEC to the BCS-side of the 
resonance. Near the resonance, the zero energy s-wave 
scattering length a exceeds the interparticle spacing and 
the interparticle interactions are unitarity limited and 
universal. Recent experiments have entered the crossover 
regime and yielded results of the interaction strength by 
measuring the cloud size and expansion. 

As in the case of bosonic clouds, the frequencies of col- 
lective modes of Fermi gases can be measured to high 
accuracy, it is of major interest to investigate their de- 
pendence on the equation of state along the crossover. 
It was pointed out yJ3 the collective frequencies of a su- 
perfluid Fermi gas at T = 0, trapped in a harmonic po- 
tential, approach well defined values in the BEC and the 
unitarity limit regimes, where the density dependence of 
the chemical potential can be inferred from general ar- 
guments. In the intermediate region, various investiga- 
tions, based on the hydrodynamic theory of superfluid 
and suitable parameterizations of the equation of state, 
have appeared recently El El El EI El El El The 
first experimental results on the collective frequencies of 
the lowest axial and radial breathing modes on ultra cold 
gases of 6 Li across the Feshbach resonance have also be- 



come available Since the BCS and the unitarity 

limits are characterized by the same collective excitation 
frequencies, there is a growing interest to study the sound 
velocity [22 . 1231 l24l | to make a clear identification of these 
two regimes and to better characterize two kinds of su- 
perfluid. 

The axial excitations of ultra cold gases in a cigar 
shaped trap can be divided into two regimes: i) long 
wavelength excitations where wavelength is equal or 
larger than the axial size, ii) short wavelength excitations 
where wavelength is much smaller than the axial size. In 
the former case, the axial excitations are discrete and 
the lowest breathing mode frequency has been measured 
[20L l2l|. In the later case, the axial excitations can be 
described by a continuous wave vector k. However, the 
finite transverse size of the system also produces a dis- 
creteness in the radial spectrum. The short wavelength 
axial phonons with different number of discrete radial 
nodes give rise to the multi branch Bogoliubov spectrum 
(MBS) H. 

The inhomogeneous density in the radial plane deter- 
mines the curvature of the mode spectrum. The effect of 
the inhomogeneous density in the radial plane decreases 
(since the radial size increases) as we go from the molec- 
ular BEC side to the weak-coupling BCS side for fixed 
number of atoms and the trapping potential. We would 
expect that the MBS will be different in the different 
regimes and it can be used to distinguish different super- 
fluid regimes along the BEC-BCS crossover. 

It should be noted that the axial excited state is cou- 
pled with the discrete radial nodes within a given angular 
momentum symmetry. For example, when we excite the 
system to study the sound propagation along the sym- 
metry axis, this perturbation inherently excites all other 
low energy transverse modes having zero angular momen- 
tum. Similarly, the above arguments are also applicable 
to other low energy mode spectrum, e. g. spectrum 
for the breathing mode. To determine the various mode 
spectrum, we must take into account that the incident of 
the mode coupling between the axial quasiparticle states 
and the transverse modes. 

In this work, we calculate the sound velocity in an 
inhomogeneous as well as homogeneous Fermi superfluid 
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along the BEC-BCS crossover. We also study the low 
energy MBS of a cigar shaped superfluid Fermi gas along 
the BEC-BCS crossover by including the mode coupling. 
It is important to study such spectrum in view of the 
current Bragg scattering experiment |2(| on the MBS of 
an elongated cloud of weakly interacting BEC. 

This paper is organized as follows. In Sec. II, we calcu- 
late the transverse eigenfrequencies and its corresponding 
eigenfunction of an elongated Fermi superfluid along the 
BEC-BCS crossover. In Sec. Ill, we discuss about the 
equation of state of the Fermi superfluid. The sound ve- 
locity, phonon mode and monopole mode spectrum are 
presented in Sec. IV. We give a brief summary and con- 
clusions in Sec. V. 



II. HYDRODYNAMIC EQUATIONS AND 
EIGEN-FREQUENCIES 

We consider a two-component Fermi gas in a long cigar 
shaped harmonic trap potential V(r,z) = (M /2){u> 2 r 2 + 
oj 2 z 2 ) at zero temperature. Here, lo z << oj r . We assume 
that the system behaves hydrodynamically throughout 
all regime. If the system is BCS superfluid, then as long 
as the oscillation frequency is below the gap frequency 
needed to break up a Cooper pair this condition is ex- 
pected to be fulfilled. The system can be described by 
the following Schrodinger equation Il6| 



iti 



9-0 



2M 



V(r) + n(n)]i>, 



(1) 



where M is the mass of the Fermi particles and /x(n) is 
the equation of state which depends on the magnitude 
and nature of the interaction strength. 

Using the Madelung transformation tp = y/ne l6 and 
neglecting the quantum pressure term, we obtain the hy- 
drodynamic equations of motion for the Fermi superfluid 
which are given by the continuity and the Euler equa- 
tions, respectively, 



and 



dn 

~dt 



M^ = -V[ M (n) 



V(r) + 2 Mv2 l 



(2) 



(3) 



Here, n(r, t) and v(r, t) = (h/M)V0 are the local density 
and superfluid velocity, respectively. We also assumed 
that u! r >> lu z so that it makes a long cigar shaped trap. 

The equation of state enters through the density- 
dependent chemical potential. We assume the power- 
law form of the equation of state as n(n) — CrO as 
in Refs. [H Q El IH IH . At equilibrium, the den- 
sity profile takes the form n (r) = {u/C) 1 / 1 (1 — f 2 ) 1 / 7 , 
where f = r/R and R = y^2u/Alu> 2 . Linearizing 
around equilibrium, n = uq + Sn, v = <5v and fi(n) = 



M( n o) + {9fi/dn)\ n — nQ 8n. The equations of motion for 
the density and velocity fluctuations are 



_ = -V-[no(r)*v], 



dSv dfi(n) 



(4) 



(5) 



Taking first-order time-derivative of Eq. (@J and using 
Eq. J^}, the second-order equation of motion for the 
density fluctuation is given by 



d 2 8n _ r , s „du(n) . . . 



8t 2 



dn 



(6) 



In the long cigar shaped trap, we assume the normal 
mode solution of the density fluctuation which can be 
written as 

Sn(r, z, t) = an(r)e i l w Wt-*«] . (7) 
Substituting Eq. Q into Eq. ©, then one can obtain 

- Cj 2 a {k)5n{r) = |v f •[(l-f 2 ) 1 /' 1 'V f (l-f 2 ) 1 - 1 / 7 5n(r) 



^~k 2 (l-r 2 )Sn(r) 7 



(8) 



where lo = uj/uj r and k = kR. Here, a is a set of two 
quantum numbers: radial quantum number, n r and the 
angular quantum number, m. 

For k = 0, it reduces to a two-dimensional eigenvalue 
problem and the solutions of it can be obtained analyti- 
cally. The energy spectrum is given by 



— \m\ + 2n r [j(n r + \m\) + 1]. 



(9) 



The corresponding normalized eigenfunction is given by 

5n a = A(l - f2)l/ 7 -l f |m|p j (l/7-l> l |)( 2 f2 _ 1 ) e im0 ! 

(10) 

where P„ ,b ' (x) is a Jacobi polynomial of order n and (f> 
is the polar angle. Also, the normalization constant A is 
given by 



A 2 _ 2 2 ~ 2/7 



[r(n r + i)] 2 r(i/ 7 )r(2/ 7 + 2n r + H) 



^R 2 r(l/ 7 - l/2)[T(l/7 + n r )} 2 T(2n r + \m\ + 1) ' 

(11) 

For 7=1, the above energy spectrum and its corre- 
sponding eigenfunction exactly matches with results of 
Ref. [27]]. Note that the modes with n r — and m =/= 
do not depend on the equation of state. This is because 
the flow in these modes are incompressible and the inter- 
nal energy does not change during the oscillation period. 
The radial breathing mode is u)\ = \/2 H + l)cj r which 
exactly matches with the result of Ref. [14j . The exper- 
imental results of the radial breathing mode pol l2*lj is 
well described ^3 this analytic spectrum. 
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The solution of Eq. JSJ) can be obtained for arbitrary 
value of k by numerical diagonalization. For k ^ 0, we 
expand the density fluctuation as 



Sn = b a Sn a (r, 



(12) 



Substituting the above expansion into Eq. JSJ, we ob- 
tain, 

= [Ql - [\m\ + 2n r (j(n r + \m\) + 1)} 



lk 2 ]b a + l~k 2 J2 M aa 'b 



Here, the matrix element M aa i is given by 



(13) 



M aa , = A 2 



x P„ 



(70, 



d 2 f{\ — f 2 ) 2 7o^2+|m| + |m'| e i(m-m')0 

m ' l \2f 2 - l) J Pi7°' |m|) (2f 2 - 1), (14) 



where 70 = 1/7 — 1. The above eigenvalue problem (Eq. 
(113(1 ) is block diagonal with no overlap between the sub- 
spaces of different angular momentum, so that the solu- 
tions to Eq. l|13l) can be obtained separately in each an- 
gular momentum subspace. We can obtain all low en- 
ergy multibranch Bogoliubov spectrum on the both sides 
of the Feshbach resonance including the unitarity limit 
from Eq. (|13|) which is our main result. Equations 1|13|) 
and (|14fl show that the spectrum depends on the average 
over the radial coordinate and the coupling between the 
axial mode and transverse modes within a given angu- 
lar momentum symmetry. Particularly, the coupling is 
important for large values of k. 



III. EQUATION OF STATE 

To calculate the sound velocity and the MBS, we need 
to know how the adiabatic index 7 depends on the two- 
body interaction strength. At zero temperature, the en- 
ergy per particle of a dilute Fermi system can be written 



e = -E F e{y), 
5 



(15) 



where Ep — h 2 k 2 F /2M is the free particle Fermi en- 
ergy and e(y) is a function of the interaction parame- 
ter y = l/kpa. In the unitarity limit (y — > O*) one 
expects that the energy per particle is proportional to 
that of a noninteracting Fermi gas. The fixed-node diffu- 
sion Monte Carlo calculation of Astrakharchik et al. psl ] 
finds e(y — ► 0) = 0.42 ± 0.01. An analogous calculation 
of Carlson et al. gave e(y -> 0)_= 0.44 ± 0.01. The 
calculation of Astrakharchik et al. |28| is quite complete 
and gives the behavior of the energy of the system across 
the unitarity limit. On the basis of the data of Carlson et 
al. [2^, Bulgac and Bertsch proposed the following 
behavior of e(y) near the unitarity limit: 



e(y) = Z-Cy- ~y 2 + o(y 3 ), 



(16) 



where £ ~ 0.44 and £ = 1 for both positive and nega- 
tive values of y. However, the data of Ref. gives 
a continuous but not differentiable behavior of e(y) near 
y = and it suggest £ = £_ = 1 in the BCS regime and 
C = C+ = 1/3 in the BEC regime. On the basis of the 
data of Astrakharchik et al. [23, Manini and Salasnich 
[l7l | proposed the following analytical fitting formula of 
e(y) for all regimes in the BEC-BCS crossover including 
the unitarity limit: 



e(2/) = oil 



a 2 tan [a 3 y 



\y\ 



ft + \y\ 



(17) 



This analytical expression is well fitted with the data of 
Ref. jl^l for a wide range of y on both sides of the reso- 
nance. We shall use Eq. (|17fl for further studies in this 
work. Two-different sets of parameters are considered in 
Ref. 0: one set in the BCS regime (y < 0) and an 
another set in the BEC regime (y > 0). In the BCS 
limit, the values of the parameters |l7j are a\ = 0.42, 
a 2 = 0.3692, a 3 = 1.044, ft = 1.4328 and ft = 0.5523. 
In the BEC limit, the values of the parameters are 
ai = 0.42, a 2 = 0.2674, a 3 = 5.04, ft = 0.1126 and 
ft = 0.4552. The advantage of a functional parameteri- 
zation of e(y) is that it allows straightforward analytical 
calculations of several physical properties. The chemical 
potential fj, is given by 



unys 

m 



/i = e(n) + 



de(n) 
dn 



E F [e{y)- y -e'{y)\, (18) 



One can extract an effective adia- 



where e'(y) - -g^ 
batic index 7 and its dependence on the scattering length 
a by defining the logarithmic derivative as |l7| 



7 = 7 



n dfj, 
li dn 



¥e'(y) 



e(y) - fe'(y) 



(19) 



The radial size of the Fermi system in all the regimes 
of the BEC-B CS crosso ver can be obtained from the re- 
lation: R = ^2u/Muo 2 r . From Eq. JTHJ, one can obtain 
the radial size which is given by 



R 



r^e{y)~ V -e'{y), (20) 



where ro = a av (24iV) 1 / 6 is the radial size of the free Fermi 
gas in a harmonic trap potential j^l), a av = yJh/Mu} av 
and w av = {us^uJz) 1 / 3 is the average oscillator frequency of 
the trap potential. In the weak-coupling BCS limit, the 
ground state energy per particle is ebcs(^) = (3/5) .Ef and 
the chemical potential is ut, cs = Ep. The corresponding 



radius is R 



bes 



,{2AN) 



1/6 



ro. In the unitarity 



limit, the ground state energy per particle is e un j(n) = 
(3/5)Ep£ and the chemical potential is u uni = Ept;. The 
corresponding radius is i? un i = a av (24A r £ 3 ) 1 / 6 = r y/^. 
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IV. SOUND VELOCITY AND MULTIBRANCH 
BOGOLIUBOV SPECTRUM 

A. Sound velocity 

Before presenting the exact numerical results, we make 
some approximation for a quantitative discussion. If we 
neglect the couplings among all other modes in the m = 
sector by setting I' = (n r ,0) in Eqs. (jT3|l and JT3J, one 
can easily get following spectrum: 

ul r = 2n r (jn r + 1) + 2(1 - M nr , nr )k 2 . (21) 

In the limit of long wavelength, the n r = mode is 
phonon-like with a sound velocity 



Ul 



(2 - 7)7 fi 



M 



(22) 



For 7 = 1, it exactly reproduces the weakly interacting 
BEC results j^. This sound velocity is different from 
the result obtained in Ref. [24|. The reason for the dif- 
ference is given below. In Ref. [24]|, the sound velocity 
is obtained by simply integrating Eq. (JHJ) on radial co- 
ordinates. In this work, we are multiplying by the com- 
plex conjugate of Sn on both sides of Eq. JSJ and then 
integrating it on radial coordinates. Since the density 
fluctuation at the lowest energy state is a function of the 
radial coordinate, the two average procedure gives two 
different result. Note that the correct average procedure 
is considered in our work. For the homogeneous Fermi 
system, the sound velocity can be obtain from Eq. 121(1 
by neglecting the M nr „ r and it is given by 



Ul 



(23) 



The sound velocity in the inhomogeneous system is 
smaller by a factor of yl — 7/2 with respect to the sound 
velocity in a homogeneous Fermi systems. This is due to 
the effect of the average over the radial variable which 
can be seen from Eqs. 1(13(1 and 1(14(1 . 

Using Eqs. fT%|> . fT§)l and (|2*2*|) . the sound velocity in 
the inhomogeneous Fermi superfluid along the BEC-BCS 
crossover including the unitarity limit is given by 



e'(y)+ v ^e"(y)}[-^+e(y)-y^e"(y)} 



My)] 



(24) 

where vf = \J2Ef/M is the Fermi velocity. Similarly, 
by using Eqs. {TSJl, (jT§)) and l(25|) . the sound velocity in 
the homogeneous Fermi superfluid along the BEC-BCS 
crossover including the unitarity limit is given by 



Ul 



n\l[\e{y)- V l e'{y) + V -e"(y)]. (25) 



Equation ((25(1 exactly agrees with the result of Ref. ^} • 



In the molecular BEC limit, the sound velocity in 
the inhomogeneous bosonic systems can be written as 
u m = y/ /i m /2M m , where \i m is the chemical potential 
of the molecular BEC and M m = 2M is the mass of 
a molecule. The chemical potential fi m can be written 
as /Lt m = 47ra m ?i 2 n„ l /M m , where n m — k F /Qir 2 is the 
molecular density and kp is the Fermi wave vector. Here, 
a rn = 0.6a is the two-body scattering length between two 
bound molecules |32| . A simple expression for the sound 
velocity in the molecular BEC limit can be written as 



V F 



0.6 1 
12ny' 



(26) 



Using equations ((24(1 . 1(25(1 and ((26(1 . we plot the sound 
velocity along the BEC-BCS crossover in Fig. 1. There is 




FIG. 1: Plots of the sound velocity along the BEC-BCS 
crossover including the unitarity limit. The solid and dashed 
lines are corresponding to the sound velocity in inhomoge- 
neous and homogeneous Fermi superfluid, respectively. The 
dot-dashed line corresponds to Eq. 12611 . 

a small kink at the unitarity limit y = due to £_ 7^ C+ • 
Otherwise, Fig. 1 shows that there is a smooth crossover 
of the sound velocity from the molecular BEC side to the 
BCS side through the unitarity limit y = 0. Fig. 1 also 
shows that Eq. 1(2611 matches very well with Eq. ((24|) fol- 
iar ge values of y. 

For homogeneous Fermi systems the sound velocity in 
the two limiting cases can be obtained from Eq. 1(25(1 
and these are given by u\ — 0.37i>f in the unitarity limit 
and by u± — 0.57vp in the weak-coupling BCS limit. 
These results exactly matches with the previous results 
[22I I23] ] . Similarly, the sound velocity for the inhomo- 
geneous Fermi systems in the two limiting cases can be 
obtain from Eq. ((24(1 and these are given by u\ — 0.30t>i? 
in the unitarity limit and u\ = 0A5vp in the dilute BCS 
limit. The sound velocity in the inhomogeneous Fermi 
system is less than that in the homogeneous Fermi sys- 
tem with the same density at the center of the trap as 
the former system. However, this difference is large in the 
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BCS side compared to the BEC side. The sound velocity 
of the inhomogeneous Fermi superfiuid can be measured 
by observing the propagation of the sound pulses along 
the symmetry axis as it is done for weakly interacting 
BEC H3- 



B. Phonon mode spectrum 

In Fig. 2 we plot the phonon mode spectrum in the 
weak-coupling BCS limit (y << 0), unitarity limit (y = 
0) and BEC side of the unitarity limit (y — 0.25) by 
solving Eq. (|13l) . These spectra have the usual form 
like bj = U\k at low momenta, where the sound velocity 
Ml is given in Eq. It is seen from Fig. 2 that 

the behavior of the phonon mode spectrum is different 
for different regimes characterizing each superfiuid phase. 
For example, the slope of the phonon spectrum in the 
BCS limit is large compared to the unitarity and BEC 
limits as expected. 



ing mode nib is 
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FIG. 2: Plots of the phonon mode spectrum in the BCS-limit 
(dot-dashed), unitarity limit (dashed) and BEC side of the 
unitarity limit with y = 0.25 (solid line). 



C. Monopole mode spectrum 

In Fig. 3, we plot the monopole mode spectrum in 
three different regimes by solving Eq. (|I3|) . In the 
long wavelength limit, the monopole mode has the free- 
particle dispersion relation with some effective mass nib 
and a gap A& = y2(7 + 1 )u>r- In the long wavelength 
limit, the breathing mode spectrum can be calculated 
from Eq. i|I 3|) by using the first-order perturbation the- 
ory. The spectrum for the monopole mode in the long 
wavelength limit is given by Wi(fc) = \/2(j + l)u> r + 
hk 2 /2nib + 0(k 4 ), where the effective mass of the breath- 



m b = M- 



r hu> r 



(2 + 7)( 7 + l) 



(2- 7 )2 



(27) 



Note that 7 = 2/3 in the BCS and unitarity limits. 
Therefore, the monopole mode frequencies are the same 
at the BCS and the unitarity limits. However, the be- 
havior of the spectrum in two different regimes are com- 
pletely different. For example, the effective mass of the 
monopole mode spectrum in the BCS limit is small com- 
pared to that of the unitarity limit. 



3.5 



2.5 - 
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k r. 







FIG. 3: Plots of the monopole mode spectrum in the BCS- 
limit (dot-dashed), unitarity limit (dashed) and BEC side of 
the unitarity limit with y — 0.25 (solid line). 



V. SUMMARY AND CONCLUSIONS 

In this work, we have calculated the sound velocity in 
the homogeneous as well as inhomogeneous Fermi super- 
fluid along the BEC-BCS crossover. The sound veloc- 
ity in the inhomogeneous Fermi superfiuid can be mea- 
sured by observing the sound pulse propagation along 
the symmetry axis, similar to the experiment by An- 
drews et al. [33 for weakly interacting BEC. The hy- 
drodynamic description presented in this work enables 
us to produce correctly all low-energy multibranch Bo- 
goliubov spectrum by including the coupling of the axial 
mode with the radial modes within the same angular mo- 
mentum sector. An analytic expression for the effective 
mass of the breathing mode spectrum is obtained. Due 
to the axial symmetry, the modes having zero angular 
momentum can be excited in the Bragg scattering ex- 
periment. Particularly, the spectrum for the phonon and 
monopole modes in the different regimes can be observed 
in the Bragg scattering experiments as these spectrum 
are observed in Ref. |2(j for weakly interacting BEC. 
By measuring the sound velocity in the pulse propaga- 
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tion experiment and by observing the low energy Bogoli- Science, 
ubov spectrum in the Bragg spectroscopy, one can make 
a clear identification of various superfluid regimes along 
the BEC-BCS crossover. 
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